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Meson correlation functions are studied in the model with four-fermion interaction Lagrangian. 

• We demonstrate that despite the singular character of system mean energy and corresponding 
quark condensate found, the meson observables are finite, quite well identified and compatible 
with experimental energy scale. It allows the similar model Hamiltonians to be used for describ- 
ing the nonequihbrium features of quark/hadron systems which reveal themselves in studying 
ultrarelativistic heavy ion collisions. The analytical results for meson correlation functions in 

■ the Keldysh model are given. 

The intensive experimental study of ultrarelativistic heavy ion collisions pushes forward a great 
interest in developing theoretical and phenomenological description of nonequihbrium processes in 

^ ', quark/hadron matter. This stage preceding a thermalization and chemical equlibration is of crucial 
importance in governing the collision process development but very complicated for reliable theo- 

^ retical interpretation. The analysis of experimental data available leads to the conclusion that at 

• the initial moment of clashing, system with great number of degrees of freedom appears and its 
constituents are strongly interacting. The characteristic time of nonequihbrium stage at the RHIC 

^ experiments is roughly estimated, for example, as 1 fm/c |lj and the energy density reached exceeds 

\ 15 GeV/fm 3 which is much higher than the corresponding quantity for the nuclear matter and its 

t— i ' value for the bag model. Clearly, such estimates urge (and allow) to speculate on the description of 
00 ' 

o 



processes at the quark level with the mechanism of dynamical mass generation included. 

In such a context the study of equilibrated states and phase diagram of quark/hadron matter is 
i-h , based on the Nambu-Jona-Lasinio model (NJL) [2] in which the adequate picture of spontaneous 
chiral symmetry breaking is properly incorporated and low energy meson physics is well understood 
[3], [1]. However, the NJL model does not accommodate the gluon degrees of freedom and is not 
directly applicable for analyzing the nonequihbrium processes. Thus, searching the related models 
which are free of such shortages and share the attractive features of NJL in the low energy region 
is a topical and practical task. In Ref. [5] the effective Hamiltonians with four-fermion interaction 
in the form of the product of two spatially separated currents mediated by a formfactor have been 
considered. As in the NJL model it is supposed that the ground state of the system is formed by the 
quark-anti-quark pairs with vacuum quantum numbers and oppositely directed momenta and in the 
framework of the Bogolyubov-Hartree-Fock approximation the description of quarks as the quasi- 
particles which is appropriate in the broad momentum range of momenta has been developed. A 
comparative analysis of the models with the formfactor behaving as the 5-function in the coordinate 
space (the NJL model) and the similar formfactor behaviour but in the momentum space (the model 
allied to the Keldysh model which is well known in the condensed matter physics) pa]) teaches that 
these dissimilar models lead to the equivalent quasi-particles when the dynamical quark masses are 
comparable. Actually, it turns out the parameters characterizing the quasiparticles are developed 
by the common dynamical mechanism which is practically insensitive to the formfactor species. 
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One unexpected feature of these models is a discontinuity of mean energy functional considered as 
a function of current quark mass what results in some difficulties at fitting the quark condensate 
beyond the chiral limit. Speaking literally the quark condensate and the mean energy of quark 
ensemble are infinite. Fortunately, neither are physically observable quantities and in order to make 
a reliable conclusion about the models we should study the meson correlation functions, for example. 
In this note we explore the meson observables in the Keldysh model which resembles a toy model 
but in some aspects elucidated below, turns out quite instructive. The meson characteristics despite 
the singular character of the ground state are finite and fairly adequate to correspond to the energy 
scale of existing experimental data. Certainly, such a conclusion does not concern all the features 
which is quite understandable within so a simple model. For example, the 7r-meson mass is slightly 
underestimated when the tuning parameters leading to the dynamical quark mass compatible with 
the NJL result are used and the pion decay constant f n disappears (see below). 



1 Model Lagrangian 

We take the model Lagrangian density (discussed in Ref. [S]) in the form of a product of two quark 
currents localized at the space coordinates x and y which are bound by the formfactor F^, i.e. 

C = q (ij^ +im) q-q t a j,, q J dy F /xv (x - y) q' t a ~f v q' , (1) 

where q = q(x,t), q = q(x,t), q' = q(y,t), q' = q(y,t) are the (anti-)quark fields, t a = X a /2 are the 
generators of SU(N C ) colour gauge group and m is a current quark mass. The Lagrangian density is 
given in the context of the Euclidean field theory and 7^ are the Hermitian Dirac matrices, — 
1, 2, 3, 4. The effective Hamiltonian corresponding to the Eq. (TjQ) results from the averaging procedure 
when the quark behaviour is strongly affected by intensive stochastic gluon field. The (anti-)instanton 
ensemble was considered as such a background. As a general form the formfactor in Eq. (pQ) can be 
presented by a sum of two components F fll/ (x — y)=G F(x — y) 5^ v + J^ u (x — y), where the second 
term is spanned on the relative distance vector. In the first component we single out the constant 
G characterizing the strength of four-fermion interaction. In particular cases when the formfactor 
F has the 5-function shape in coordinate space we come to the NJL model. With the formfactor 
behaving as F(p) = (2tt) 3 S(p) the model is similar to the Keldysh model. In order to simplify the 
consideration we ignore the contribution of the correlator corresponding to J^ij))- The unprejudiced 
analysis of the system behaviour beyond the chiral limit performed in Ref. [5] demonstrates that the 
four-fermion interaction develops a singularity. The mean energy of the ensemble goes to infinity 
and the quark condensate demonstrates the singular behaviour as well. In the leading order in the 
iV c -expansion we obtain for the generators of colour group J2 a =i ^tj^li ~ 2 "^"fcj anc ^ utilizing the 
Fierz transformation 7 M (8> 7/x = 1 <8> 1 + 275 <8> 275 ~ §7/x <8> 7/x ~ §7^75 ® 7^75' we nave f° r the scalar 
contribution in the mean field approximation the following effective Lagrangian density 

£ ~ q (il^ + im) q- J dy G F(x - y ) (q q') q' q . (2) 

where the angle brackets denote the corresponding average. It is interesting to notice the interaction 
term is composed with the colourless quark operators with x and y coordinates interchanged. The 
selfconsistency condition allows us to extract the dynamical quark mass as 

M(p) = 2N c [-%GF(p- q ) m + M{q) m . (3) 

At the spontaneous breaking of chiral symmetry takes place at this stage we have for the formfactor 
behaving as F(x) = 5(x) the well known gap equation 

M = 2N r G / W ^% m + M 1/2 , 

(27T) 3 [q2 + (m + M)2] l/2 
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where A NJL is the cut off parameter. For the Keldysh model in the mean field approximation it looks 
like 

M(p) = 2N C G m + M(p) 1/2 . (4) 

Transforming this solution into the function p(M) we have for the quark mass in the chiral limit 
in the Keldysh model M(p) = [(2N C G) 2 — p 2 ] 1 ^ 2 - In the further analysis bearing in mind the Fierz 
transformation we consider a more general form of the Lagrangian with two constants included. 
They characterize the interacting strength for the scalar and pseudoscalar channels — G and for the 
vector and the axial-vector ones — G v . Clearly, the scalar and pseudo-scalar coupling constats are 
identical in the chiral limit [2]. 

2 Bosonization 

Here we introduce the meson fields adapting well known bosonization procedure for the scalar channel 
(similar relations are valid for the other channels). It is convenient to introduce the auxiliary variables 

Q s (x,y;t) = F 1 ' 2 {x-y) - [ q{x;t)q{y;t) + q{y;t)q{x;t) } , 

1 (5) 
Q a ( Xl y;t) = F l ' 2 {x-y) - [ q(x;t)q(y;t) - q(y;t)q(x;t) } , 

for the symmetric and anti-symmetric combinations of quarks. Since the formfactor is a symmetric 
function with respect to an interchange of coordinates x — > y then interaction contribution in the 
scalar sector which is just the point of our interest can be written as 

G 1 G 

vL = Y F ( x - y) 2 t *) *) + ^ *) t)q{y; t) ) = - {q 2 s - qD . 

It is easy to see that now the standard procedure of bosonization may be realized with the Gaussian 
integration which concerns auxiliary meson fields a s (x,y;t) and a a (x,y;t). Indeed, the integration 

a fG\ 1/2 

is performed over the combinations including the meson and quark fields as . — ( ^ ) Qs, 



(2cy 

G\ 1/2 

' • Q a . Then the interaction term may be presented in the form including the meson 



l<Jg _ 

fields as 

vL = j(Q 2 s - Ql) - 

- + F" 2 {x - y) q(x; t)q(y; t) + F^{x - y) q(y; t)q(x; t) . 

and integrating over the quark fields we obtain the effective theory operating with the mesons only. 
We do not show the detailed calculations here because they follow the standard procedure with 
one minor distinction which is a doubling of meson fields. Because of the same reason we present 
the succinct exposition of calculating the equation for dynamical quark mass and extracting meson 
correlator behaviour. We remind only that the first variation of the effective action allows us to 
determine the dynamical quark mass 

-i- + <Q,) = 0, -? T + l (Q a ) = 0, (7) 
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(here the first equation for symmetric field is equivalent to the Eq. (JHJ)). As a result for the induced 
quark mass we havaS 

M(x) = F^(x) <jM{x) , a?>{x - y) = -~ F^{x - y) (<j\q(x)q(y) + q{y)q{x)\a) . (8) 

Certainly, we are interested in the real solutions and should consider only the case when the average 
contribution of anti-symmetric quark combination becomes trivial (Q a ) = 0. 

For the quadratic terms of effective meson Lagrangian in the scalar channel we have 



2G 



+ \ q(x; t) F^(x - y) q(y; t) + q(y; t) F^(x - y) q(x; tj) ■ 

^ q(x'; F'/\x' - y>) q(y'; t>) + q(y'; t>) F^(x> - y>) q(x'; t>) 



here we imply that the meson fields with the primes are dependent on the coordinates x' and y' and 
t' . Then the pairing of quark fields with utilizing the corresponding Green functions leads in the 
momentum representation (the integrations over the corresponding 'internal' variables are dropped) 
to the equation 

0-aO, Q, Pi) K a a p(p, q, p', q', p 4 ) ap(p', q' , p 4 ) = (9) 
^(p,9;p 4 K(-p -9;-P4) + ^(p, g ;p 4 K(-p -95-P4) 1 TlS{ k; h )S(k - p - q;k 4 - Pa ) ■ 



2G 2 



0"'o - «cr' 



2 v r ' v 2 

+ (Ts ~ icra F 1 ^(k - q)F 1 ' 2 {k + <z')^4^ + as ~ 2 l ° a F l / 2 {k - ^^(fe + p' f' s ~ 2 i<T ' a 

where the following notations are introduced a a , s = cr a , s (p, Q', Pi), ^' as = a a , s (p', Q 1 ', — Pa)] q' = ~P — 
q — p' , a, j3 = s, a. 

Apparently, there is no special need to investigate the meson correlation functions in so general 
form. We can obtain quite enough information on the solutions analyzing some particular cases. First 
we consider the conditions when the formfactors becomes identical F 1 ^ 2 {k — q) = F 1 l 2 {k—p — q—p'), 
F 1 ^ 2 (k — q) = F l l 2 {k + p'). These allow to conclude that the momenta of quarks coincide p = q, 
i.e. there is no a relative motion of quarks in such a situation. Then it is easy to understand that 
the contribution of antisymmetric fields becomes degenerate in this configuration and the remaining 
symmetric part corresponds explicitly to the standard bosonization procedure. Finally, we have for 
the meson correlators in scalar and pseudoscalar channels 

K a,, = _± +2Nc [ F {k-p) hih ~ Pi) + ~ 2P) T{m + M ^ m + M ( fc - 



2G C J (2vr) 4 v ^' [k 2 + (m + M(k)) 2 }[(k i - Pi ) 2 + (k-2p) 2 + (m + M(k-2p)) 2 y 
k 2 = k\ + k 2 . In particular, for the Keldysh model it reads 

r,a,n _ 1 | o at [ dk ± fc * ( fc4 ~ ~P 2 T(m + M{p)f 

2G +ll " c J 2n [k 2 + E 2 (p)][(h-p,) 2 + E 2 (p)] ' [W) 

where the notation E 2 (p) = p 2 + (m + M(p)) 2 for the quark energy is used. To simplify the 
presentation of formulae we omit the energy E dependence on the momentum. Due to the fact that 



1 In the momentum representation, respectively, we have M(p) = G j , ^ 3 F(p — q) i Tr S(q), where S(q) is the 

(27Tj 

quark Green function. 
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only the integration over k± is essential we say about the one- dimensional model for mesons in this 
paper. 

The denominator of the 7r-meson in Eq. (II 01) can be written in more convenient form as 



[hi + E 2 } {(h - Pa? + E 2 ] [2 h (k 4 -p A ) +P l + 2 E 2 



1 



+ 



1 



k 2 + E 2 (k 4 -p 4 ) 2 + E 2 
We see the correlation function of the 7r-meson in Eq. ( TTUl) is expressed by three integrals 



^ = __ +2JVc/ *i 



+ 



1 



2tt 2 \k 2 + E 2 (k 4 -p 4 ) 2 + E 2 
Calculating them we have finally the following result 

pI + Ap 2 



dkd 



pl + Ap 2 



2tt [k 2 + E 2 ][(k 4 -p 4 ) 2 + E 2 )] 



Kn = -JG + JE +N < < 



2p 4 E (p 4 + 2iE) ' 
1 _ p\ + Ap 2 
2E E {p{ + AE 2 ) ' 
p\ + Ap 2 



Ira p 4 > iE , 
| Im p 4 \ < E , 
Im p 4 < —iE 



2 p 4 E (p 4 - 2iE) ' 

where p = \p\. In the Euclidean domain we have for the real values of energy p 4 that 



' E 



E (p 2 + AE 2 ) \G 



PZ + 



E-G 



G 



AE 2 + Ap 2 



(12) 



where G = 2N C G. Then we find that the meson correlation function resembles a screening factor. 
In order to investigate the pseudo-euclidean situation we continue the p 4 variable to the imaginary 
axis. Introducing the notation p 4 = iPq we have 



K* = Nr 



G 



E (P 2 - AE 2 ) 
P E (A + 2E) \ a 



G 

E d2 , G-2E 



P 2 4- 



G 



E P + 2p< 



Pn < E, 



Pn > E, 



Comparing this expression at Pq < E with Eq. ( fl2l) we make certain that the transition from the 
Euclidean variables p 4 to pseudo-euclidean ones do not change its form. The branch Im p 4 < —iE 
is not be considered because of the symmetry reason. 

Now we continue with searching the 7r-meson dispersion law which is defined by the zeros of 
correlation function K 71 = 0. The results for scalar and pseudo- scalar mesons are presented in Fig. 
1 and for the vector and axial vector mesons are shown in Fig. 2. For the branch Pq < E the 
dispersion can be received from the following equation 



It can be obtained from Eq. 



P 2 = A(E-G)E + ^(2p) 2 . 

for the quark energy 

m + M 



E = G 



M 



(13) 



(14) 



In particular, for zero quark momentum we have M(0) = G for induced quark mass and for the 
7r-meson energy we receive 

P 2 = A (E - G) E = A m (m + G) , (15) 
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which means that in the chiral limit (m — > 0) the Goldstone theorem is valid (see. Fig. 1, curve 1). 
Tuning the model parameters as was proposed in Ref. [5| we get the following parameters for the NJL 
model [3J m = (m M +m^)/2 = 5. MeV, G = 286 MeV. For the 7r-meson energy we have at that Pq « 76 
MeV for the zero quark momentum but the quark energy at low momenta is E « m + M ~ 286 



(E-N C G) 2 + ^U 



MeV. When P > P the interesting branch looks like P = N C G — E + 

but the analysis shows this branch does not satisfy the constraint Pq > E. 

Turning now to the scalar channel we present the integrand in the convenient form as (TTUT) 



1/2 



k 4 {k 4 — p 4 ) — P 2 



[(kj + E 2 ]ih- Pi y + E< 



k 2 + E 2 



+ 



1 



+ P 2 



1 



p\ + 4P 2 



2[kj + E 2 ][(k i -p 4 ) 2 + E< 



(16) 



and for the a-meson correlation function receive 



1 



+ 2N C 



dki 1 



2G ' ~""J 2vr 2\k 2 + E 2 (fc 4 - p 4 ) 2 + E 2 / 
Calculating the integrals we come to the following result 



iV, 



dfc 4 pj + 4P 2 

27 [A; 2 + P 2 ][(A; 4 -p 4 ) 2 + P 2 



K° 



1 , iV c , /v 
~2U + 2P +7Vc 



p 4 - 2iP 
' 2p A E 
1 

~2P ' 
p 4 + 2iE 

2 p 4 E 



Im p 4 > iE, 
\Im p 4 | < P, 
im p 4 < — zP. 



In the pseudo-euclidean regime for the branch P < P the a-meson correlation function is degenerated 



and for Pn > E we have P' 7 



— 2f3p — ' Comparing to the quark energy of Eq. ( 1T4"1) 



2T7 ^lt-^o 
we conclude that at P > there are not the interesting zeros in K a . 

We have considered the configuration when the relative momentum of quark and anti-quark equals 

zero and below we address the quark and anti-quark system with zero total momentum p + q = 0, 

see. Eq. (Q. For the outgoing quark momenta two configurations are possible: a) p' — p and b) 

p' = —p. In the a)-situation we obtain for the correlation functions in scalar and pseudo-scalar 

channels 



n * Ki p Tip 



1 vr + 7T 



2N r 



2 2G ' " c 2 
2 2G + c 2 



dk 4 k 4 (k 4 - p 4 ) +p 2 + (m + M) 2 
~2t7 [k 2 + E 2 }[(h-p,) 2 + E 2 ] 
dk 4 k 4 (k 4 - Pa) + p 2 - (m + M) 2 
~2j7 [k 2 + E 2 }[{k 4 - p 4 ) 2 + E 2 ] 



where a, (3 — s, a. Then for the 7r-meson we have 



-Po + G- 2E 
2G(P + 2E) 
P 2 - IE {E - G) 
2G(AE 2 



Po>P 



Pu<E 



-Po — G — 2E 



Pn>E 



2G(P + 2E) ' 
P 2 - 4P (P + G) p p 
2G{AE 2 -P 2 ) ' n< * 



The dispersion law for the Kg correlation function at Po < E is extracted from the following equation 

4(P-G)P, (17) 



P 2 



which is in a full agreement with Eq. ( 1131) if the total momentum of quark anti-quark pair is taken 
to develop value 2p — > (see. Fig. 1, curve 2). For the K£ correlation function there is no solution 
meeting the constraint Pq < E and there is no any solution for the PJ correlation function at Po > E 
as well as for the branch PJ. 
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In the scalar channel one can obtain 




-EP 2 + EP (G - 2E) + 2G{m + M) 2 [ -EP 2 - EP (G + 2E) - 2G[m + Mf 

2GEP ^P + 2E) ' = I 2GEP^P Q + 2E) 

EP 2 - AE 2 (E — G) -4G(m + M) 2 a ) EP 2 - AE 2 (E - G) + 4G(m + Mf 
2GE 2 (AE 2 - Pj) ' { 2GE 2 {AE 2 - PJ) ' 



The upper ratios are written for P Q > E and the lower ones for P < E. The dispersion law for the 
branch K s a at Pq < E is determined by the solution of the following equation 

P 2 = A {E-G) E + 4 % {m + Mf , (18) 

E 

see the dashed curve in Fig. 1. We did not manage to find the appropriate solution of dispersion 
equation for the correlation function K° in this case as the condition Pq < E is invalid. Analysis 
of the correlation function roots at P$ > E gives the same message that the suitable solutions are 
absent. We omit the discussion of the configuration b) because the results already given demonstrate 
how rich and complicated the analysis of solution branches could be. We would like to mention 
only that the presence of bound state even for the quarks with comparatively large momenta looks 
improbable. 

500 ii i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 



400 - 




p MeV 



Figure 1: The cr-meson (dashed curve) and 7r-meson (curves 1 and 2) energies in MeV as the functions 
of momenta. The curve 1 corresponds to the configuration of zero relative quark momentum. The 
curve 2 describes the quark and anti-quark system with zero total momentum. The dots demonstrate 
the quark energy E. 

In order to calculate the pion decay constant it is necessary to calculate a loop integral which is 
similar to Eq. (jTOl) in which one of the vertices responsible for the weak interaction of the quarks 
does not contain the formfactor F l l 2 {k — p) relevant for the meson fields. Then for the Keldysh 
model the integral with weak singularity develops formally the zero value what leads to the pion 
decay constant equal to zero. 

One can consider the vector and axial-vector meson channels in the similar way if the correspond- 
ing substitutions qq — > q^^q, a — > V^, qi^q — > qj^^q, vr — > are done in the relevant formulae. 
However, we omit those calculations here and do use the method of the equations for vertex functions 
(Bethe-Salpeter equation) to analyze the correlation functions quantatively. 
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3 Vertex functions 

It was demonstrated above that the Gaussian integration should be performed with the symmetric 
and anti-symmetric combinations of auxiliary meson fields and the anti-symmetric fields should 
contain an imaginary unit factor. Clearly, the corresponding analysis of the functional integral 
saddle point in the imaginary space appears to be a rather complicated task. Hence, it would be 
highly desirable to consider some alternative possibilities of studying the meson observables keeping 
under control the quark degrees of freedom. Here we are going to use the method based on the 
Bethe-Salpeter equation. We start calculating an effective quark interaction and sum up the set of 
some special diagrams as follows 





+ ••• 

It should be taken into account from the beginning that each line of these graphs has to be depicted 
as the doubled one because we consider the nonlocal meson fields. The first diagram of this set 
describes the initial interaction T {p,q) = (2n) 3 G T F{p — q) where the matrix T denotes the 
interaction channel 1, 275, 7 M , 757^. If we single out the combination q a {x;t) T a/ 3{x,y;t) qp{y;t) 
then the following equation can be calculated for the series sum 



-(2tt) 3 G T F{p -q) + T(-p, -q, r 4 ) = 

r dk A dk ~ 
: - / — G T F(-k - q) Tr S(k; fc 4 ) T S{k - p + q; k 4 - r 4 ) T(k-p + q, k; h ~ u), 



where r 4 = p 4 — g 4 and p 4 , g 4 are the corresponding components of quark and anti-quark momenta. 

We search for the solution for the vertex function in the Keldysh model, for example, in the form 

qr;r 4 ) = (2tt) 3 F 5(p — q) V(p;r 4 ) assuming that for the imaginary values of r 4 the solution 

for the vertex function possesses the pole singularity. Picking out the singular contributions we can 

obtain approximately that — G + V = G U V where the polarization operator in the Keldysh model 

dk 

can be represented in the form n = / Tr S(—p; fc 4 ) T S(—p; k^ — r^) T. Therefore we have for the 

G 

vertex function V = ^ -q and its denominator zeros determine the pole positions. Discussing 
the bosonization above we have already calculated the similar polarization operators. Omitting the 
intermediate calculaions we show here the results for pseudo-scalar and scalar channels at |r 4 | < E 
as 

1 - G n ps _ AE (E - G) - Rl 

G ~ G {AE 2 - Rj) ' 
1 - G n s _ {AE 2 -R 2 ) E + A[{m + M) 2 - E 2 ] G 



G G E {AE 2 - R t 



where R = P — Q , P = ip±, Q = iq^. In the pseudo-scalar channel the 7r-meson dispersion law 
coincides explicitly with Eq. ffTTj) and for the cx-meson we have Eq. ( 118ft . 
For vector and axial-vector channels in the Keldysh model we have 



n v, A ^_ 47V f dh h{k 4 - r 4 ) - p 2 =f (m + M) 2 
1144 ' " 4iVc J 2tt \k\ + E 2 } [(h - n) 2 + E 2 } ' 

nY: A -4iv {2h - n) Vl (20) 



dki k^ik^ — r 4 ) dij — 2 piPj + Sij p 2 ± Sij {m + M)' 
2^ [kl + E 2 ][{h-r 4 ) 2 + E 2 ] 
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II 4 - = Yl i4 . It is easy to see that these results coincide with the corresponding meson correlation 
function obtained above. Calculating the integrals in vector and axial- vector channels for Im r 4 > %E 
we have: 



1J -44 — z,iv c ? 1J -44 



2N C i r 4 E - 2p 2 



r 4 ' E r 4 (r 4 + 2iE) ' 



= -2N C ^- , n A = nj , (21) 

= i r 4 ff % - 2p, Pj = [i r 4 E - 2(m + M) 2 ]<% - 2p iPj - 

ij c r 4 {r 4 + 2iE)E ' ij ' E r 4 (r 4 + 2iE) 

At |r 4 | < E we obtain 



TT V — D TT A 



8iV c {m + Mf 



44 ' 44 £ r 2 + 4£ 2 ' 

ni = o, 1^ = 14, (22) 

n v _ 9 N E 2 Sij-PiPj tta _ 8N c P 2 5ij - PiPj 
v c E (r 2 + AE 2 ) ' ij E r\ + AE 2 ' 

Now we diagonalize the correlation functions using the fact that corresponding quadratic forms 
determine simply the Lagrangian of free vector and axial- vector mesons 

iT-A = Qv, A y2 + 2 £v, A ^ + ^. £v, A ^ ^ (23) 

where we imply the summation over the indices which are repeated. It is valid by definition that 
= $iu> ~ n^ A . If we redefine the space components of (axial-) vector fields by substituting 

Vi = Vi + a PiV 4 , where a v ' A = ~ ^ % , and exclude the mixed components V 4 Vi from quadratic 

Pi ( \i Pj 

form ( BSD we § e ^ that the fourth component of vector field enters the quadratic form with coefficient 



(C Y - A n ) 2 

C 44 A = CJ 4 A — - — 4 ~ v ^ . The components of the tensor Cj- A remain unchanged. The numerical 

Pi C'ij Pj 

analysis demonstrates the acceptable solution for the dispersion of the fourth component exists 
only for the axial-vector field at Rq > E. However, we do not discuss this solution in the present 
paper. The spatial components of vector fields are searched as having two different forms for the 

transversal Vi = (^5ij — vj- and longitudinal Vi = Pi t>" components. First we consider the 

case of Rq > E. The dispersion law for the transversal component of vector field v 1 - has the form 
Rq = G Y — 2E, where G v = 2N C G V . Therefore the nontrivial solution is possible when the condition 
G v > 3E is satisfied. For clarity we take the constant as G v = 1.5G which corresponds at the low 
quark momentum to the value G v ~ 1.5E and therefore for the fitting parameter set selected this 
solution branch does not manifest itself. The dispersion of longitudinal component v" in this case is 

(E - N C G Y ) 2 + 4iVc g vP . The branches of 

the vector and axial-vector channels are depicted in Fig. El The branch discussed above is shown 
by the solid curve and denoted by 1. For the transversal component of axial- vector meson a 1 - we 

\ E - N cGv ) 2 + 4^G v (m + M)' 

dashed line and marked by 2. The dispersion of the longitudinal component a" is calculated from 
Ro = G v . This branch appears when the condition G v > E is satisfied. In Fig. [2] it corresponds 
to the dashed straight line 3. We see at low momentum the longitudinal a" and transversal a 1 - 
components practically coincides. 

Now we address the situation of Rq < E. The dispersion of the transversal component of the 
vector meson v x is defined by -Rq = AE (E — G v ). This curve is shown in Fig. El by solid line 4. 



defined as R Q = N C G V — E + 



obtain R = N C G V -E + 



21 1/2 



. This curve is depicted in Fig. El by 
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Figure 2: The relative energies Rq of vector (solid curves) and axial- vector (dashed curves) mesons 
(in MeV), as the functions of momenta (in MeV). The notations are explained in the text and the 
dots denote the quark energy E. 



For the longitudinal components v " we have Rl = 4 [iP-G ^Y j, 

and the solid curve 5 

shows its behaviour. For the transversal component of the axial-vector meson a 1 - we obtain Rq = 

AE 2 — 8N C G v ^r. The corresponding curve is presented as the dashed line 6. Its right-hand 
component is practically degenerated with the curve 4 because the induced quark mass goes to zero 
at large momenta and the restoration of chiral symmetry takes place. The longitudinal component 
of axial- vector field a" is degenerated. 

The 7r-meson energy in the NJL model for the tuning parameter values considered coincides with 
experimental data and looks like E^ 3h = 140 MeV. In Ref. [5] it was supposed that for the Keldysh 
model the dynamical quark mass in the low momentum region is equal to the dynamical quark mass 
of the NJL model. This assumption has led to the almost identical quasiparticles for both models. 
But the 7T- meson energy turned out rather underestimated E^ = 76 MeV although the relative scale 
of mesons for different channels was maintained. The axial-vector meson had gotten heavier than 
the vector meson and the meson of pseudoscalar channel was the lightest one, indeed. 

Calculating the correlation functions here we did not exploit the convenient trick of shifting the 
integration variable to make the integrals symmetric and were keeping the integration contour fixed. 
The one-dimensional character of the model provides us with the obvious possibility to study the 
dependence of the correlation functions on the integration contour shape as well. In particular, it 
is interesting to trace the turn of integration contour to the imaginary axis. Then two quark poles 
could also be treated as the "Wigner" phase and the continuation of correlation functions will be 
different from what has been done in this note. 



4 Correlation functions in the Minkowski space 

The model developed allows us to study easily the meson correlation functions in the Minkowski 
space as well and to compare them to what we obtained above. In fact, the task is technically 
related to computing the following integral (dependent on the Euclidean variables) within the fixed 
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contour 



f°° dq 
J-oo 2ix 
r°° dq 
J-oo 2% 



(q 2 + E 2 ) [(q — p) 2 + E 2 



dq 1 
2~ix (2iEf 



2iE 



1 



+ 



2iE 



1 



+ 



p(p — 2iE) q — iE p(p + 2iE) q — p — iE 
-2iE 1 -2iE 1 



+ 



p(p + 2iE) q + iE p(p — 2iE) q — p + iE 



which is the sum of four pole terms: iE designated as ei, p + iE indicated as e 2 , — iE marked as m\ 
and p — iE signed as m 2 . It is implied the parameter p may take the complex values as well. Then 
the calculation of this integral leads to 

l!t In) p > //•; . 



—I 
AE 1 



p(p + 2iE) ' 

AiE 
p 2 + AE 2 ' 
2iE 
{ p(p - 2iE) ' 



\Im p\ < E , 
Im p < —iE 



Our concern here is the particular situation when the parameter p is pure imaginary p = iP. Now we 
should calculate the similar integral with another fixed contour which corresponds to the Minkowski 
space (turned to 90 degrees regarding the Euclidean integration) contour 



J 



f°° dq 

J-oo 2^ (q 2 - 
00 dq i 
-oo 2^ [2{E - is) 



(q 2 -E 2 + ie){(q-p) 2 -E 2 + ie] 



2(E-ie) 



+ 



2(E - is) 



p[p - 2(E - is)] q- E + ie p[p + 2(E - is)] q - p - E + is 
-2(E-ie) 1 -2(E-ie) 1 



+ 



+ 



p[p + 2(E - ie)\ q + E -ie p[p - 2{E - is)] q-p + E 



IE 



and its calculation gives 



J 



1 

AW 



2E 

p\p-2{E-ie)\ 
AE 

p 2 - A(E - ie 
2E 

{ p[p + 2(E-is) 



71 5 



Im p > ie , 
\Im p\ < e , 
Im p < —ie . 



Compared to the Euclidean configuration the pole contributions to the integral in the Minkowski 
space are interchanged, i.e. the poles of e-type become the poles of m-type. Thus, the result ob- 
tained for the Minkowski space will be valid for the Euclidean configuration if the contribution of 
the m 2 -pole at Im p > iE is omitted and the contribution of e 2 -pole at Im p < —iE is added. One 
should not change anything at \Im p\ < E. The figure shows one of possble integration contours 
corresponding to the situation Im p > iE. Similarly at Im p < —iE the contour should be deformed 
in order to have the contribution of e 2 -pole. The detailed analysis make possible to formulate the 
general rule for reproducing the proper result in the Minkowski space. The integration contour should 
be deformed in such a way to have the contributions of the e± and e 2 - poles only, i.e. the integration 
contour looks like being squeezed in between the poles e 2 and m 2 . 
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= e 2 



ei 



rn 2 



mi 



In order to calculate the integrals / and J in ^-representation we rewrite the integral / (using 
the well-known identity and implying the analytical continuation of all auxiliary functions in the 
parameter p) as 



dq 



dx 



dq 



oo 2tt (q 2 + E 2 )[(q - pf + E 2 } Jo J-oo 2tt {x(q 2 + E 2 ) + (1 - x)[(q - pf + E 2 }} 2 



Ap 3 L 



1/2 
-1/2 



E 2 
p 2 " 



3/2 



E(p 2 + 4E 2 ) ' 



where ^ — x — 1/2. Such an integral treatment makes transparent that the result above corresponds 
to the calculation of / in the fixed contour for the branch when \Im p\ < E. Changing the variables 
of integration as E 2 — > — E 2 we are able to reproduce the corresponding result for J. Thus, we may 
conclude that calculating in the ^-representation fully rereproduces the result for the Minkowski 
space being analytically continued. 

Considering the meson correlation functions in the Minkowski space we are interested in the 
situation when the external parameters are real what corresponds to the constraint \Im p\ < e for 
the J integral. It is a pretty simple task to obtain the final results for the corresponding dispersion 
laws using the results of calculations for the Euclidean space. Below we show as an example the 
results for the (anti- ) quark total momentum equal zero 

P 2 = AE(E-G) , P 2 =AE 2 - A^p 2 , 
P 2 ± = 4E(E - G v ) , P 2 A± = AE 2 - 4^ p 2 , 
P 2 V] =AE 2 -A^ (m + M) 2 , 



the axial vector field correlator a\\ becomes degenerate. At G — G v 7r-meson becomes degenerate 
with the vector meson and cr-meson with the axial vector meson. It is clear the realistic relations 
between the meson masses correspond to the situation when G v < G. In order to give another 
example we take G v = G/2 (in addition to the Euclidean consideration in which the constant was 
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Figure 3: The meson energies obtained by using the correlation functions in the Minkowski space 
as the functions of quark momentum. The double quark energy is depicted by dots. The dispersion 
law for the longitudinal component of vector meson field is depicted by stars and the dashed line 
with the symbol v± corresponds to the transverse component. The dashed line with the symbol a± 
shows the dispersion law for the transverse component of axial vector meson. The solid line with the 
symbol a is devoted for the scalar meson whereas the similar line with the symbol ir corresponds to 
7r-meson. 

taken as G v = 1.5 G. The meson energies as quark momentum functions are shown in Fig. [3] in 
the Minkowski space. Curiously, the bound states of qurk and anti-quark do exist at any quark 
momentum for the present configuration. 

5 Conclusion 

In this note we demonstrate that despite the singular behaviour of mean energy and quark condensate 
which was observed in Ref. [5] the meson observables are finite, well identified and compatible with 
the experimental energy scale. The number of effective degrees of freedom which define the quasi- 
particle picture in the NJL and Keldysh model are comparable. The Keldysh model being as simple 
as the NJL model looks like a suitable candidate for describing the nonequilibrium processes in the 
(anti-) quark ensembles. Due to the one- dimensional character of the Keldysh model the analytical 
continuation from the Euclidean region of meson observables to the pseudo-euclidean one is easily 
performed and controled. The amazing feature of our consideration is that the bound states are 
revealed at any quark momenta in Minkowsky space. Eventually we conclude that if the quasi- 
particles in the different models are similar the meson observables are also alike. 
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